Effect of serum sample storage temperature on metabolomic and proteomic biomarkers

Prospective biomarker studies can be used to identify biomarkers predictive of disease onset. However, if serum biomarkers are measured years after their collection, the storage conditions might affect analyte concentrations. Few data exists concerning which metabolites and proteins are affected by storage at − 20 °C vs − 80 °C. Our objectives were to document analytes affected by storage of serum samples at − 20 °C vs − 80 °C, and to identify those indicative of the storage temperature. We utilized liquid chromatography tandem mass spectrometry and Luminex to quantify 300 analytes from serum samples of 16 Finnish individuals with type 1 diabetes, with split-aliquot samples stored at − 80 °C and − 20 °C for a median of 4.2 years. Results were validated in 315 Finnish and 916 Scottish individuals with type 1 diabetes, stored at − 20 °C and at − 80 °C, respectively. After quality control, we analysed 193 metabolites and proteins of which 120 were apparently unaffected and 15 clearly susceptible to storage at − 20 °C vs − 80 °C. Further, we identified serum glutamate/glutamine ratio greater than 0.20 as a biomarker of storage at − 20 °C vs − 80 °C. The results provide a catalogue of analytes unaffected and affected by storage at − 20 °C vs − 80 °C and biomarkers indicative of sub-optimal storage.


Scientific Reports
| (2022) 12:4571 | https://doi.org/10.1038/s41598-022-08429-0 www.nature.com/scientificreports/ freeze-thaw cycles. This might inadvertently lead to including low-quality samples, or, excluding valuable samples in biomarker discovery and validation studies due to excessive caution, resulting in failure to identify potential biomarkers or limiting the use of a particular bio-resource. Many pre-analytical factors, i.e., those processes before sample measurement, potentially affect metabolomic and proteomic studies in serum samples. These include study design, sample collection, sample handling and storage, and sample preparation 3,4 . For example, up to 3% of the detected metabolome was affected by multiple freeze-thaw cycles or/and extended thawing at 4 °C 5 . Surprisingly little data are found on the effect of temperature during long-term storage of serum samples, and more specifically, on the effect of storage at − 20 °C vs − 80 °C. To the best of our knowledge only studies investigating short-term storage are available 6,7 .
The objectives of this study were two-fold: in serum samples collected from subjects with type 1 diabetes, first, to identify biomarkers of storage at − 20 °C vs − 80 °C indicative of sub-optimal storage in general; and second, to define serum metabolites and proteins unaffected by storage at − 20 °C vs − 80 °C.

Materials and methods
Participants and serum samples. The Finnish Diabetic Nephropathy (FinnDiane) Study is a prospective nationwide multicenter study comprising more than 8400 adults with type 1 diabetes from 21 university and central hospitals, 33 district hospitals, and 26 primary health care centers across Finland covering more than 10% of all the individuals with type 1 diabetes in Finland 8 . Although, strictly speaking, it is not a populationbased study, the geographical distribution of the FinnDiane patients closely follows the distribution of the general population of Finland. Patients participated in the study during a regular visit to their attending physician during which detailed demographic and medical history data were collected with standardized questionnaires.
The Scottish Diabetes Research Network Type 1 Bioresource (SDRNT1BIO) 9 is a prospective cohort study, comprising 6127 people with a clinical diagnosis of type 1 diabetes, representative of all adults with type 1 diabetes in Scotland, and recruited between December 2010 and November 2013. At recruitment, clinical measurements and a blood sample were taken.
Both studies were performed in accordance with the Declaration of Helsinki; all participants gave their written informed consent, and the study protocol was approved by the local ethics committees ( For this study, we selected a sub-population of 16 FinnDiane individuals, for whom split-aliquot serum samples, stored continuously at − 20 °C and − 80 °C and never thawed, were available. This storage dataset was used as the discovery dataset. The study was performed as a part of a larger biomarker project which included 315 individuals from the FinnDiane and 916 individuals from the SDRNT1BIO cohorts 10 . The individuals in the larger biomarker project were used to validate biomarkers of storage at − 20 °C vs − 80 °C, as the serum samples from FinnDiane were stored at − 20 °C and the serum samples SDRNT1BIO were stored at − 80 °C. In SDRNT1BIO it took on average 2 h 29 min between when the sample was taken to when it entered the freezers after processing (median time 2 h 15 min, interquartile range: 1 h 35 min-3 h 13 min) and they were not thawed prior to analysis. We do not have data on time from sample collection to freezer for FinnDiane and most of the 315 samples were thawed at least once before analysis.

Biomarkers.
Altogether 300 analytes were measured in non-fasting serum samples using two different targeted platforms: 269 (122 metabolite concentrations and metabolite ratios, and 147 tryptic peptides (relative quantitation)) were measured by liquid chromatography tandem mass spectrometry (LC-MSMS) at the Well-Child Laboratory (King's College London, UK) and 31 protein concentrations were measured using the Luminex platform at Myriad RBM (Austin, TX, USA) (Supplementary Table 1).
Several quality control steps were applied to the data before the main analysis. If an analyte had only a single constant value (n = 28) or was missing more than 50% in the 32 measured discovery samples from 16 individuals (n = 8), it was removed from the analysis as uninformative. In addition, 2 analytes failed quantification: ApoD (518.3/824.8) (precursor ion m/z/product ion m/z) due to variable chromatographic separation, while the Tissue Inhibitor of Metalloproteinases 3 (TIMP-3) assay did not meet analytical specifications.
We evaluated the reproducibility of the analyte quantification by measuring duplicate samples (35 duplicate samples in SDRNT1BIO, 25 in FinnDiane) in the pilot phase of the main biomarker project. We calculated intra-class correlation on the duplicate samples in both cohorts separately and removed analytes with intra-class correlation < 0.4 (n = 79) in either cohort.
After removing 107 analytes from the analysis, 193 analytes (27 proteins, 83 metabolites and metabolite ratios, and 83 tryptic peptides) were included. The analyte median (IQR) values or the removal reason in the discovery dataset are listed in Supplementary Table 1.
Left-censored values (Supplementary Table 1 and Ref. 10 ) below the detection threshold were imputed to half of the detection threshold and right-censored values (not present in the discovery samples and only present for N-terminal prohormone of brain natriuretic peptide in the validation samples) were imputed to maximum reported value.

Statistical analysis.
To quantify the effect of storage at − 20 °C on the serum analytes we defined a score Given that samples were paired and within each pair the only difference between samples was the storage temperature, a paired t-test or a paired Wilcoxon signed rank test was performed to quantify the statistical significance of the pairwise difference. The paired Wilcoxon signed rank test was used if the distribution of the pairwise differences did not seem to follow a normal distribution (Shapiro-Wilk test p-value < 0.05), otherwise the paired t-test was utilized.
Further, to determine the effect of storage time on the stability of the analytes a linear model was fitted: where fc bm is the ratio of analyte levels at − 80 °C and − 20 °C (fold change) and t stor is the storage time in years. Fold change was used to quantify the effect of storage temperature on the analyte levels to obtain a measure independent of the mean magnitude of the measurement pair. Through the log 2 transformation, we are modelling the effect of storage time expressed in terms doubling (or halving) of the ratio between analyte levels at the two storage temperatures.
To identify biomarkers predictive of storage temperature, the individual analyte's ability to discriminate between − 20 °C and − 80 °C samples was evaluated by calculating the receiver operating characteristics (ROC) curve and the corresponding area under the curve (AUC) in the discovery set, using the raw analyte values. In other words, only the raw analyte value was used as the explanatory variable to predict if a sample was stored at − 20 °C or − 80 °C and a ROC curve was constructed using this simple classifier for each potential biomarker. The discriminative power of biomarkers classifying − 20 °C and − 80 °C samples correctly in the discovery set was tested in the combined FinnDiane-SDRNT1BIO validation set. Further, the discovery set was used to find for each biomarker the interval that separates the − 20 °C and − 80 °C samples, such that setting the classification threshold within this interval gives perfect classification in the discovery set. Then ranges of sensitivity and specificity were calculated in the validation set when the classification threshold was within this interval.
Biomarker panels. We also constructed parsimonious biomarker panels to predict the sample storage temperature using LASSO-penalized regression implemented in the glmnet R-package 11 , where the optimal value of the penalty parameter was learned through internal leave-one-out cross-validation. We first constructed panels separately for the proteins, metabolites and tryptic peptides using the FinnDiane discovery dataset. Based on the historically well-known instability of glutamine in stored serum samples 12 , we also constructed a panel for the metabolites in which glutamate/glutamine ratio was forced into the model. The performance of all panels was tested in the two validation sets by applying the regression coefficients learned on the discovery dataset.
Missing analyte values were imputed to median and analyte values were log 10 -transformed prior to constructing the biomarker panels.

Analytes affected by storage at − 20 °C vs − 80 °C.
Out of 193 analytes passing quality control (27 proteins, 83 metabolites and metabolite ratios, and 83 tryptic peptides) there were 12 analytes affected by − 20 °C vs. − 80 °C storage temperature, defined as | S d | > 1.5 (p < 2.6 × 10 -4 for each analyte, corresponding to a Bonferroni corrected significance threshold of p = 0.05 for multiple testing of 193 analytes). Eight analytes had S d > 1.5 , indicating markedly higher concentration when stored at − 80 °C, namely kallistatin (643.4/971.6), neutrophil gelatinase-associated lipocalin, methionine, free cystine, glutamine, C4 carnitine, C2 carnitine and sulphocysteine. Conversely there were four analytes for which S d < −1.5 (p < 2.6 × 10 -4 ) indicating lower concentration when stored at − 80 °C, namely glutamate, free carnitine, interleukin-1 receptor type 2 and interleukin-1 receptor type 1 ( Table 2). On the contrary, 120 analytes showed only minimal effect of storage temperature with a |S d | < 0.5 and were defined as robust (Supplementary Table 2). The distribution of S d across all analytes is shown in Fig. 1. www.nature.com/scientificreports/ For eight analytes the storage time was associated with log 2 (fc bm ) at a Bonferroni corrected significance threshold of p = 0.05/193 = 2.6 × 10 -4 . In other words, the ratio between the paired − 80 °C and − 20 °C samples changed over time for these analytes: C2 carnitine, C4 carnitine, C3 carnitine, sulphocysteine, glutamate/glutamine, free cystine, C6 carnitine and glutamine ( Table 2 and Supplementary Fig. 1). Glutamate/glutamine ratio had the biggest change with an estimated fold change of 0.67 per year; for the seven other analytes, the concentrations increased by a fold change ranging from 1.09 to 1.41 per year. To illustrate, for these seven analytes with fold change > 1, the ratio between the paired − 80 °C and − 20 °C samples was higher in the samples that were stored for a longer time compared to the ones stored for a shorter time.
Biomarkers of storage at − 20 °C vs − 80 °C. The biomarkers discriminating perfectly between samples stored at − 20 °C and − 80 °C in the discovery dataset were free cystine, glutamate, glutamate/glutamine ratio and kallistatin (643.4/971.6). These biomarkers also separated the samples in the validation dataset almost perfectly (AUC > 0.997); glutamate/glutamine ratio had an AUC of 1 ( Table 3). The distributions of the log 2 -transformed biomarker values differed also in visual inspection (Fig. 2). The intervals separating the − 20 °C and − 80 °C samples in the discovery dataset are given in Table 4 Biomarker panels. The four biomarkers for storage at − 20 °C vs. − 80 °C were detected with LC-MSMS metabolite and tryptic peptide platforms. As they may not be readily available in different biomarker studies, we further investigated biomarkers for the storage temperature separately on each biomarker platform. We identified parsimonious biomarker panels separately for proteins, metabolites and tryptic peptides using LASSOpenalized regression (Supplementary Table 3). The panels for each platform consisted of one to five biomarkers and performed well in classifying the samples stored at − 80 °C from the samples stored at − 20 °C: AUC > 0.988 for all panels in the discovery dataset and AUC > 0.995 for all the panels in the validation datasets (Table 5 and  Supplementary Table 4).

Discussion
Metabolomic and proteomic analyses of patient cohort samples are central to biomarker discovery studies aimed at developing new clinical diagnostics and prognostics. It is generally considered that biomarker discovery is relatively straightforward for acute clinical presentations, as they are usually accompanied by significant protein and/or enzyme release, a characteristic that has been used to develop many of our current repertoire of clinical diagnostics. In addition, rapid outcome measures allow informative samples to be obtained within a relatively short time window and, therefore, should not require long term storage before analysis.
Less straightforward, but arguably clinically more valuable, are the recent initiatives to elucidate new biomarkers that can be used to either define the risk of developing a chronic disease or enable very early detection of evolving disease processes. The aims are to ensure that, in high-risk individuals, early therapeutic intervention ameliorates disease progression, and, in low risk individuals, reassurance is afforded, and clinical intervention minimised. In addition, in the current era, it is anticipated that new clinical biomarkers will aid cost-effective assessment of new therapeutics. Unfortunately, the progression of chronic diseases, e.g. cirrhosis, tends to be slow and the development of clinical complications insidious, e.g. diabetic nephropathy. The implication is that, to identify early diagnostic/prognostic biomarkers in slowly progressing diseases requires the collection and storage of biological samples many years or even decades in advance of the manifestation of easily measurable clinical endpoints. This is particularly true regarding the clinical complications of both type 1 diabetes and type 2 diabetes. In the case of type 1 diabetes the diagnosis is usually made very early in the disease, providing the opportunity to collect biological samples, to establish the natural history of the development of disease complications, both cross-sectionally and longitudinally. If these samples are going to be valuable in future metabolomic and proteomic biomarker studies, it is essential that they are stored appropriately.
On the contrary, many of the studies that collected samples a decade or more ago, stored them at − 20 °C, and the samples may have been subject to multiple freeze-thaw cycles and other sub-optimal pre-analytical factors and storage conditions. To assess the usability of such studies in the current biomarker analyses, we evaluated  Table 2). The www.nature.com/scientificreports/ implication is that, it would be possible, where samples from a highly informative clinical cohort have been suboptimally stored, to use the cohort samples in discovery and validation studies, provided the metabolites and/ or proteins have been previously demonstrated to be stable. Furthermore, we defined a serum glutamate/glutamine ratio > 0.20 as being highly sensitive and specific in identifying samples that have been stored at − 20 °C compared to − 80 °C. While we only considered the storage temperature in this analysis (a pre-analytical factor that is often known), we hypothesize that the glutamate/ glutamine ratio can be used as a general indicator of sub-optimal storage conditions due to other, unknown pre-analytical factors.
In recent years considerable effort has been expended in defining the effects of sample collection tubes, initial processing before prolonged storage, shipment, storage conditions, and freeze-thaw cycles on metabolomic 12,13 and proteomic profiles 13 . Studies on the effects of long-term sample storage have focused on − 80 °C 14,15 . However, there were no internationally recognised sample collection and storage guidelines in place when many of the major clinical cohort studies were originally instigated and, primarily due to cost considerations, many of the samples were stored at − 20 °C. In the case of plasma/serum samples stored continuously at − 20 °C the temperature is insufficient to ensure complete "freezing", resulting in slow but continued enzymatic conversion of metabolites and protease breakdown of proteins [16][17][18][19] . In addition, there are many non-enzymatic processes that continue in samples frozen at − 20 °C, e.g. Schiff base formation, dehydration, hydrolysis, carbamylation, and oxidation 20 , all of which will not only alter specific metabolite concentrations 21 but also, almost certainly, modify the immunospecificity of certain proteins 22 . Freeze-thaw cycles are even more critical as they provide significant opportunities for further enzymatic and non-enzymatic metabolism. Finally, even in screw-topped tubes there will be significant sublimation of samples, if they are stored for a decade or more at − 20 °C.
Although not investigated in this study, it is appropriate to point out that the temperature of storage of urine samples is even more critical, particularly for proteins: significant losses of all proteins occur in a high proportion of urine samples during initial freezing at − 20 °C 23 .
A consideration that most biological matrix metabolite and protein stability studies fail to consider is the role of the underlying clinical condition. The majority of stability studies are based on relatively few samples from healthy subjects and the results are then applied to clinical cohort samples on the false assumption that all samples are equal. Plasma enzyme activities vary significantly depending on the disease process being investigated and the underlying sample matrix. In diseases with known extreme phenotypes, e.g. renal disease and uraemia, it is critical that metabolite and protein stability are investigated using appropriate clinical samples. Diabetes is the most obvious extreme phenotype with, at least, the effects of both hyperglycaemia/glycosuria, ease of Schiff base formation with reactive glucose 21 , and increased lysosomal enzyme release 24 , enzymatic conversion and modification, to be considered.
The reasonable concerns over biological sample stability has resulted in significant improvements in current best practice guidelines for optimal collection and storage of plasma/serum samples, including splitting samples into multiple aliquots for single use, to avoid multiple freeze-thaw cycles, and storage at − 80 °C 2 . The LacaScore has been suggested as a quality control of optimal pre-centrifugation time and temperature of plasma samples 25 . In a similar fashion it would be valuable to define an objective quality control check of optimal prolonged sample storage. Table 4. Ranges of candidate biomarkers in the discovery set samples stored at − 20 °C and − 80 °C. Classification interval is the interval which gives perfect classification results in the discovery set. Sensitivity and specificity ranges are given in the validation set when the classification threshold is in the perfect classification interval. www.nature.com/scientificreports/ Although we had only a limited number of paired samples to compare storage data at − 20 °C and, consensus best practice, storage at − 80 °C, the differences between certain metabolites were so dramatic that it is relatively easy to establish a criterion for sub-optimal storage. It is interesting to note, though not unexpected, that the serum concentrations of the three sulphur-containing amino acids, free cystine, methionine, and sulphocysteine, are significantly lower when stored at − 20 °C. This is likely to be primarily due to oxidation, though in the case of free cystine the formation of protein linked disulphides could be a major contributory factor 26,27 . The reduced short chain acylcarnitine concentrations are the result of hydrolysis and are, consequently, mirrored by an increase in free carnitine. However, it is the decrease in glutamine and parallel increase in glutamate concentrations that are most informative in determining sub-optimal storage. Glutamine is considered to be relatively stable in aqueous solutions 28 but at higher temperatures, where losses can be considerable, there is no equivalent increase in glutamate. In this instance it is likely that the product, due to cyclisation, is pyroglutamate 29 . In contrast, it is well recognised within the clinical plasma/serum amino acid analysis community that optimal measurement of glutamine requires fairly quick separation of serum from whole blood and immediate analysis or storage of the serum at − 80 °C 30 . In addition, the decrease in glutamine is usually associated with a significant increase in glutamate, suggesting that, at least in part, the reaction is enzymatic, with glutaminase being the most likely candidate. The current data are consistent with conversion of glutamine to glutamate and pyroglutamate in serum samples stored at − 20 °C (Supplementary Table 2): the sums of glutamate, glutamine, and pyroglutamate concentrations are equivalent, approximately 600 µmol/L, in the samples under both storage conditions.
The associations between glutamine and glutamate/glutamine ratio with storage time at − 20 °C suggest that this ratio may be used to define varying degrees of sub-optimal storage. In this small dataset, inspection of the serum glutamate/glutamine ratio and simply taking the mean of the highest ratio in the samples stored at − 80 °C and the lowest ratio in those stored at − 20 °C provides a pragmatic discriminatory ratio of 0.20. When applied in the validation dataset the AUC was 1, demonstrating the sensitivity and specificity of this ratio in providing an objective and practical method of determining sub-optimal storage of serum samples. Although this ratio has only been validated in serum samples from individuals with type 1 diabetes, it is highly likely, given the history and our understanding of the conversion process, that it will be applicable to serum samples from disease cohorts other than diabetes. Consequently, we would argue that the glutamate/glutamine ratio should be measured in all serum samples used in metabolomic and proteomic studies as a quality control biomarker of sub-optimal sample storage.
Free cystine and the tryptic peptide of kallistatin (m/z 634.4/971.6) also provide virtually perfect discrimination and, where available, could support the glutamate/glutamine ratio.
Serum retains significant protease and peptidase activities, so it might be expected that storage of samples at − 20 °C would result in significant losses of a number of proteins and tryptic peptides of the 27 proteins and 83 tryptic peptides studied. Interestingly, this was not the case, with only kallistatin (643.4/971.6) and neutrophil gelatinase associated lipocalin (NGAL) particularly susceptible to loss of concentration. In contrast, 14 (52%) proteins and 69 (83%) tryptic peptides were robust against storage temperature. An oblique confirmation of this observation is provided by the metabolite stability data, as any significant proteolytic activity in stored samples would result in increased free amino acid concentrations; this was not the case. Difficult to explain is the apparent increase in the 2 interleukin receptors though this may be due to enzyme activity releasing or exposing more receptor.
To conclude, this study indicates just how vulnerable metabolomic and proteomic cohort studies are to sub-optimal sample storage. We defined a serum glutamate/glutamine ratio > 0.20 as being highly sensitive and specific in identifying samples, from individuals with type 1 diabetes, that have been stored at − 20 °C compared to − 80 °C, and suggest that the ratio can be used as an indicator of sub-optimal storage conditions in general. In addition, we reported 15 serum metabolites and proteins that are highly susceptible to storage at − 20 °C vs − 80 °C and provide an evidence base for the exclusion from untargeted biomarker discovery studies of serum samples stored at − 20 °C. On the contrary, we also identified 120 of the 193 examined serum metabolites and proteins relatively unaffected by sub-optimal storage, supporting that also cohort samples stored at − 20 °C may be used in biomarker studies.

Data availability
Summary level data is available in the supplementary information. Individual level data cannot be shared due to patient privacy.